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Abstract 

In continuation of our systematic effort to understand hadronic matter at high density, we 
study dense skyrmion matter and its chiral phase structure in an effective field theory imple- 
mented with the trace anomaly of QCD applicable in the large Nc limit. By incorporating 
a dilaton field x associated with broken conformal symmetry of QCD into the simplest form 
of skyrmion Lagrangian, we simulate the effect of "sliding vacua" influenced by the presence 
of matter and obtain what could correspond to the "intrinsic dependence" on the background 
of the system, i.e., matter density or temperature, that results when a generic chiral effective 
field theory of strong interactions is matched to QCD at a matching scale near the chiral scale 
47r/7r ~ 1 GeV. The properties of the Goldstone pions and the dilaton scalar near the 
chiral phase transition are studied by looking at the pertinent excitations of given quantum 
numbers on top of a skyrmion matter and their behavior in the vicinity of the phase transition 
from Goldstone mode to Wigner mode characterized by the changeover from the FCC crystal to 
the half-skyrmion CC crystal. We recover from the model certain features that are connected to 
Brown-Rho scaling and that suggest how to give a precise meaning to the latter in the framework 
of an effective field theory that is matched to QCD . 

Pacs: 12.39-x, 13.60.Hb, 14.65-q, 14.70Dj 
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1 Introduction 



In trying to understand what happens to hadrons under extreme conditions, e.g., at high density 
as in compact stars or at high temperature as in relativistic heavy-ion colhsions, it is necessary 
that the theory or model adopted for the description be consistent with the basic tenet of 
QCD. What this means in terms of effective theories using macroscopic degrees of freedom, i.e., 
hadrons, the effective field theory should be matched to QCD at a scale close to the chiral scale 
A-^ ~ 47r/7r ~ 1 GeV. In a recent important development reviewed in pP, Harada and Yamawaki 
show how this matching can be effectuated in the framework of hidden local symmetry (HLS) 
theory. In this theory, the light-quark vector meson fields figure as gauge fields coupled to 
the Goldstone pion fields vr. When this theory is matched to QCD via current correlators at a 
suitable scale Km ~ A-^, among the variety of flow paths by which the theory can follow to its 
multiple fixed points as the scale is changed, it picks the particular flow that leads to one unique 
fixed point as one reaches the point at which the phase change to Wigner mode from Goldstone 
mode takes place. This fixed point corresponds to the "vector manifestation (VM)" in which 
the parameters of the Lagrangian have the limiting behavior: 

5^0, a^l (1) 

where the parametric g is the (hidden) gauge coupling constant and a is the ratio Ft^/F^ where 
i^TT and Ffj are respectively the parametric pion ^ and would-be Goldstone scalar decay decay 
constants. In Harada and Yamawaki discuss how the fixed point is reached when one 
dials the number of flavors to a critical value A^^^; > 3 at which the physical pion decay constant 
/tt vanishes. Subsequently it has been shown that chiral phase transition at critical temperature 
Tc [HI and at critical density ric |3 involves the same VM. The consequence is that at the critical 
point, independently of whether it is driven by temperature, density or large N^, the gauge 
boson mass (both parametric and pole) approaches zero and so does the constituent quark 
mass. Thus all light-quark hadrons (except for the Goldstone pions) are expected to become 
massless in the chiral limit in some power of the gauge coupling constant g. This observation led 
Brown and Rho to conjecture that the property of VM and Brown-Rho scaling are intimately 
connected [Hj- An important lesson one can learn from this development is that without the 
Wilsonian matching to QCD, there is no way to pick the right one from the multitude of different 
ways to describe hadron properties near chiral restoration. This means that the theories that 
possess all relevant symmetries but are unmatched to QCD have no predictive power for hadron 
properties. This caveat applies not only to HLS theory but also to all chiral field theories with 
or without vector-meson degrees of freedom. 

The matching at Km of the correlators between the HLS sector and the QCD sector gives the 
parameters of the "bare" HLS Lagrangian in terms of the quantities figuring in the QCD sector, 
namely, the color gauge coupling a^, Kqcd and quark and gluon condensates. In medium, 
the condensates depend on the background, that is, on density ^. Imposing that the quark 
condensate {qq) ~ as an order parameter - vanish at the critical density ric, one obtains, at 
the matching scale, the conditions on the HLS parameters g{KM,nc) = and a{KM,ric) = 1. 
However it gives no constraint on the parametric F.jj{Km, ric)- Now bringing the scale from Km to 
the scale appropriate to the physics involved is done by renormalization group equations. Thus 
the properties both at the matching scale in matter-free vacuum and at the VM in medium 
are completely determined by the theory. But the HLS theory formalism does not offer any 

^We will denote the physical pion decay constant by the lower case to be distinguished from the parametric 
one. 

^Similar arguments hold in heat bath although we do not specifically mention it in what follows. 
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method to describe how the parameters behave at an arbitrary scale A and density n away 
from the two special points. This means that we know practically nothing on the intrinsic 
density dependence of the parameters of the Lagrangian and hence cannot compute the density 
dependence of physical quantities except at n ~ and n ~ ric. Thus it would be futile to 
attempt to employ HLS/VM theory in its present form for analyzing experimental data which 
of course sample all ranges of density. 

In this paper, we tackle the intrinsic dependence and related issues with the help of the unified 
approach to high density developed in [7]. In [7j, it was shown using a skyrmion Lagrangian 
in its simplest form that an effective Lagrangian valid in the large Nc limit (where Nc is the 
number of colors) can describe both infinite nucleonic matter and pionic fiuctuations thereon 
as density is increased and that by describing the chiral restoration as a changeover from an 
FCC crystal to the half-skyrmion matter, one can describe the properties of various excitations 
commensurate with the background given by the half-skyrmion crystal configuration. In this 
approach, the density dependence of the background is taken into account to all orders. No low- 
density approximation whose validity is in doubt except at very low density is ever made in the 
calculation. The power of the approach is that the dynamics of the background and excitations 
thereon can be treated in a unified way on the same footing with a single Lagrangian. 

In addressing the problem at hand, the Lagrangian used in [2j is probably incomplete. In 
fact, it is not clear that the intrinsic density dependence required by the matching to QCD 
discussed above is fully implemented in the model. One puzzling feature we found in |2j was 
that the Wigner phase represented by the half-skyrmion matter with (trC/o) = supported a 
non-vanishing pion decay constant. This was interpreted there as a possible signal for a pseudo- 
gap phase. However it can also be interpreted as an analog to Georgi's "vector limit" in which 
chiral symmetry is restored with pions present with non- vanishing pion decay constant. Harada 
and Yamawaki argue ^ that this phase is not consistent with QCD Ward identity, implying 
that it cannot be realized in nature. 

Our proposal in this paper is that we can circumvent the above difficulty and achieve our 
objective within a large Nc framework by incorporating into the skyrmion description the trace 
anomaly of QCD which is lacking in the standard skyrmion model. In fact it is this feature that 
led to BR scaling in j2j when the problem was treated albeit in a schematic way. In this paper 
we formulate BR scaling of ^ in a more rigorous way and offer a scheme that could shed light 
on the intrinsic density dependence needed in the HLS/VM theory. 

The basic idea in our approach is as follows. We consider a skyrmion- type Lagrangian with 
spontaneously broken chiral symmetry and scale symmetry associated, respectively, with nearly 
massless quarks and trace anomaly of QCD. Such a theory may be considered as an A^c — > oo 
approximation to QCD. Suppose that the Lagrangian can describe not only the lowest-excitation, 
i.e., pionic, sector but also the baryonic sector and massive vector meson sector all lying below 
the chiral scale A^. The skyrmion description gives not only the single baryon spectra but also 
multi-baryon systems including infinite nuclear matter. We are interested in how low-energy 
degrees of freedom in many-body systems behave in dense matter, in particular as the density 
reaches a density at which QCD predicts a phase transition from the broken chiral symmetry 
to the unbroken chiral symmetry or chiral restoration in short. For this purpose, one first 
looks for the ground state of the many-baryon system in question as a soliton solution of the 
Lagrangian and then looks at the fluctuation of effective fields in various channels of low-energy 
excitations with the effect of the background taken into account self-consistently. As one varies 
the density of the system, the parameters of the theory involved in the process must adapt to 
the density of the skyrmion background. The problem we are interested in is whether we can 
extract from the unified scheme an information on the intrinsic dependence inherent in effective 
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field theories matclied to QCD , i.e., HLS/VM theory of Harada and Yamawaki. The hope is 
that by looking at the structure at the phase transition and tuning the parameters such that 
the features predicted by HLS/VM are reproduced, we can ultimately learn about the intrinsic 
dependence of the parameters and hence "derive" BR scaling. We will see that some progress 
in this direction can be made. 

The content of this paper is as follows. In Section 2, we give our model Lagrangian which is 
the simplest form of skyrmion Lagrangian, namely, the original Skyrme Lagrangian implemented 
with trace anomaly of QCD. In Section 3, a single skyrmion is analyzed to define the parameters 
of the theory at zero density. The chiral phase transition from an FCC crystal to a half-skyrmion 
CC crystal is described in terms of the model Lagrangian in Section 4. Section 5 describes 
how the pseudo-Goldstone pion and the scalar of the trace anomaly ("dilaton") behave in the 
skyrmion matter as a function of density. In Section 6, the chiral restoration transition is made 
from the inhomogeneous phase - to which the FCC crystal collapses - to the half-skyrmion 
crystal phase. The contact with BR scaling, and indirectly with the parametric property of 
HLS/VM, is made in this section. Some concluding remarks are given in Section 7. 

2 The Model Lagrangian 

The starting point of our work is the skyrmion Lagrangian introduced by Ellis and Lanik [S] 
and employed by Brown and Rho [5] for nuclear physics that incorporates the trace anomaly of 
QCD. Here we will study the same Lagrangian from a more modern perspective. 

The classical QCD action of scale dimension 4 in the chiral limit is scale-invariant under the 
scale transformation 

x ^ = X'^x, A > 0, (2) 

under which the quark field and the gluon fields transform with the scale dimension 3/2 and 
1, respectively. The quark mass term of scale dimension 3 breaks this scale invariance. At the 
quantum level, scale invariance is also broken by dimensional transmutation even for massless 
quarks, as signaled by a non- vanishing trace of the energy- momentum tensor. Equivalently, this 
phenomenon can be formulated by the non- vanishing divergence of the dilatation current D^, 
the so called trace anomaly, 

d^D^ = = E ^1^9 - ^1VG'^.G"^^ (3) 

q 9 

where l3{g) is the beta function of QCD. 

We will implement broken scale invariance into large Nc physics by modifying the skyrmion 
Lagrangian, 

C = ^Trid^U^d'^U) + -^Tr[C/ta^C/, U^d.U? + ^^Tt{U + [/t - 2). (4) 
4 32e^ 4 

The chiral field U = exp(ir • vt/Ztt) has scale dimension and therefore the Lagrangian respects 
neither the scale invariance of QCD nor its breaking eq. (jSJ- The current algebra term with two 
derivatives is of dimesion 2 and the meson mass term is of dimension 0, while the Skyrme term 
with four derivatives has scale dimension 4. 

In order to make the Skyrme model well behaved under the scaling properties of QCD , we 
introduce an additional degree of freedom in the form of a scalar field x with a scale dimension 
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1, whose coupling to the U fields is defined by, 



X^(ln(x//x)-^)+^ 
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We have denoted the non vanishing vacuum expectation value of x as J-^, a constant which 
describes, as we shall see, the decay of the scalar into pions. The second term of the trace 
anomaly (jSJ can be reproduced by the potential energy V{'x) for the scalar field, which is 
adjusted in the Lagrangian © so that V = dV/dx = and d^V/dx^ = m\aix = fx- 

The vacuum state of the Lagrangian at zero baryon number density is defined by ?7 = 1 and 
X = fx- ^il^ take the latter to be positive and therefore the field values around the vacuum 
will be positive. The fiuctuations of the pion and the scalar fields about this vacuum, defined 
through 

U = exp(if • (f/ /^), and x = fx + X (6) 

give physical meaning to the model parameters: /jr as the pion decay constant, as the pion 
mass, fy. as the scalar decay constant, and as the scalar mass. For the pions, we use their 
empirical values as 7,^ = 93MeV and m.,^ = 140MeV. We fix the Skyrme parameter e to 4.75 
from the axial- vector coupling constant g^i as in ref. jlUj . However, for the scalar field Xi 
experimental values for the corresponding parameters are available yet. 

The scalar field may be interpreted as a bound state of gluons, the so-called glueball ^ . If one 
assumes that it is a pure glueball, then one can use the gluon condensate Go = {0\{as/iT)G^uG'^'^\0) 
to restrict the product of its mass rriy. by the vacuum expectation value f^ as 

l/x"^x = (|f?o)i (7) 

Fluctuating around the vacuum at zero density, the Lagrangian yields the relevant interactions 
for the process x ~^ In the chiral limit, the coupling from the current algebra term is given 
by 

~ 3 

Cx.n = fY.i9M'- (8) 

a=l 

This yields the decay width r(x ^ vrvr) as 

3m?, mt 1 

in the x rest frame, where we used eq. (jTj) to eliminate in favor of Gq. Introducing the 
pion mass into the calculation modifies the expression slightly (8j. Using the ITEP value Gq = 
0.012GeV^ [n]^ one obtains 

r(x ^ vrvr) ^ 0.6GeV x (m^/lGeYf. (10) 

This equation provides a restriction to the range of the possible gluonium masses 

r(x ^ vrvr) < 1 MeV for < 0.3 GeV 
r(x ~^ Trvr) > 0.5m^ for > 1 GeV. 



^Tho more recent ITEP value [HI is Go = 0.011 ± 0.009 GeV" which is consistent with the older value. 
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Table 1: Parameter sets of the model Lagrangian 







fx 


Go 


set A 


550MeV 


240MeV 


0.004 GeV^ 


set B 


720MeV 


240MeV 


0.007 GeV^ 


set C 


lOOOMeV 


240MeV 


0.012 GeV^ 



We see that a light scalar gluonium with the mass below 400 MeV cannot be detected because 
of its too narrow decay width, while a heavy one with the mass greater than 1 GeV cannot be 
recognized as a well-defined resonance. 

In the following development, we will find that the scalar x ™- this picture cannot describe 
pure gluonium; it must contain a light scalar quarkonium component. This "soft" component 
cannot vanish in the chiral broken phase in order to satisfy the symmetry properties of the 
fundamental theory ^Sj- The low-energy sum rules|n|, related to the gluon condensates, are 
controlled by the gluonium component and can be used to fit the parameters of the theory in 
the free case, i.e. and from Gq. There can however be a subtle "mended symmetry" that 
puts the mass of x nearly degenerate with that of the p meson, i.e., ~ 700 MeV |12j . It is not 
clear what happens in dense medium. Our conjecture JS| is that as the density is increased, the 
two components "decouple" with only the quarkonium component interacting with matter. We 
can think of this as the "hard" gluonium component being integrated out in an effective field 
theoretic sense with their properties absorbed into the definition of the fields and the parameters 
jlfij . This makes sense as the characteristic mass scale drops compared with the mass scale of 
the gluonium. Chiral symmetry restoration implies the vanishing of the "soft" component. 

There are some indications that the above conjecture makes sense. In ref. ^HI; the scalar 
field is incorporated into a relativistic hadronic model for nuclear matter not only to account for 
the anomalous scaling behavior but also to provide the mid-range nucleon-nucleon attraction. 
Then, the parameters and rriy^ are adjusted so that the model fits finite nuclei. One of the 
parameter sets is = 550 MeV and = 240 MeV (Set A). On the other hand. Song et 
al. |19j obtain the "best" values for the parameters of the effective chiral Lagrangian with the 
"soft" scalar fields so that the results are consistent with the "Brown- Rho" scaling ; explicitly, 
niy. = 720 MeV and = 240 MeV (Set B). These observations are consistent with the notion 
that "mended symmetry" is operative for the scalar field to be identified with the sigma field of 
linear sigma model in certain "dilaton" limit 

For completeness, we consider also a parameter set of rriy. = 1 GeV and fy. = 240 MeV (Set 
C) corresponding to a mass scale comparable to that of chiral symmetry A,^ ~ ^T^f-w- 

The parameter sets defined above for the model Lagrangian are summarized in Table 1. 

3 Single Skyrmion with the Scalar Field 

The topological baryon number current associated with the homotopy of the mapping UQ{f) : 
53(i?3 _ {oo}) ^ ^3 of gjj^2) is defined as 

i?^ = ^e''''^f'Tv{U'od,UoU^odxUoUldpUo). (12) 
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The soliton solution with the baryon number B = 1 can be found by generalizing the spherical 
hedgehog Ansatz of the original Skyrme model as 



Uo{r) = exp(zf ■ fF(r)), and xo{r) = fxC{r), (13) 



with two radial functions F{r) and C(r). 

Then, the mass of the single soliton can be expressed as 



Mb=i = 47r / r'^dr 



flr^2(j,,2 , ^sin^FX 1 siv? F ( sir? F 2 



i^C'[F" + 2^^\+—^^^[^^ + 2F 



^flmlC\X - cosF) + ^-^ ic'^ + ^{(C^lnC - i) + 1} 



.(14) 



where the prime means the derivation with respect to r. Variations of Mb=i with respect to 
F{r) and C(r) lead to the following equations of motion 



fyC^ + '-^)F"ir) + 2fA rC^ + )F'(.) - 



2 /2C2 sin F . cos F - ^'^^^^"^""'^ - f^lC^' sin F = 0, (15) 



for F(r) and 



f2c"ir) + ^C'(r) - /2C(r) (^F'\r) + ^^^] - 3/>2 (1 _ ^osF)C\r) 

e2 2ri3 



+/>^C^(r)lnC(r) =0, (16) 
for C(r). 

At infinity, the fields Uo{f) and Xo(^) should reach their vacuum values. The asymptotic 
behavior of the equations of motion reveals that the radial functions reach the corresponding 
vacuum values as 

F{x) — , and C{x) ~ 1 — . (17) 

In order for the solution to carry a baryon number, Uq has the value —1 at the origin, that is, 
F(x = 0) = vr, while there is no such topological constraint for C{x = 0). All that is required is 
that it be a positive number below 1. The equations of motion tell us that for small x, 



F{x) ~ TT - ax + + O(x^), 
C{x) ~ Co + fix"^ + 0{x 



2 , r^t^^^ (18) 



The coupled equations of motion for F{x) and C{x) together with the boundary conditions 
H17|) and (|18|) can be solved in various different ways. For example, we can use an iteration 
method. We first start with Cj=o(x) = 1 for all x and solve the equation of motion (|15|) with 
C{x) fixed as Cq{x) to obtain Fj=o(^)) which is nothing but the profile function of the single 
skyrmion of the original Lagrangian. Then, we take this Fo(x) as F{x) in the equation of motion 
(|TH|) to obtain Ci=i[x). We repeat this iteration until Fi+i(x) and Ci+i(x) converge to Fj(x) 
and Ci{x). This iteration converges quite fast and after ten iterations the solution can be found 
with a sufficient accuracy. 



6 



1.00 



0.75 



C(X) 



0.50 



set A: MB.i=1364MeV, <t^>^'^=o.44fm 
:1389MeV, <t^>^'==0.44fm 



set B : M 
set C : M 



B=1 



= 1416MeV,<t^>''^=0.43fm 



0.25 



0.00 I ■ ■ ■ ■ ' ■ ■ ■ ■ ' ■ ■ ■ ■ ' ■ — ■ ■ ■ — I 

0.0 0.5 1.0 1.5 2.0 

r(fm) 

Figure 1: Profile functions F{x) and C{x) as a function of x. 

Shown in Fig. ^ are the profile functions F{x) and C{x) as a function of x(= efj^r). F{r) 
and consequently the root mean square radius of the baryon charge 

{r^/^ = (^J dVB^{f)Y^ (19) 

show little dependence on m^. On the other hand, the changes in C(r) and the soliton mass 
are recognizable. The larger the scalar mass is, the smaller its coupling to the pionic field 
and the less its effect on the single skyrmion. In the limit of rriy. — > oo, the scalar field is 
completely decoupled from the pions and the model returns back to the original one, where 
C(r) = 1, Msoi = 1479 MeV and (r^)^/^ = 0.43 fm. The strong dependence of C(r) on 
may come from the asymptotic behavior (|17j) at large x. As for the solion mass, note that Mgoi 
scales approximately as {f-K/e) and C^(r) is multiplied to in the current algebra term of the 
Lagrangian. Thus, C(r) < 1 reduces the effective fj^ inside the single skyrmion so that the 
soliton mass decreases accordingly. For example, for the parameter set A, the soliton mass gets 
7% reduction and we can imagine that the effective /tt is reduced by the same amount in average 
in the region where the single skyrmion is located. 

4 FCC Skyrmion Crystal with the Scalar Field 

Consider a crystal configuration made up of skyrmions, where each FCC lattice site is occupied 
by a single skyrmion center with Uq = —1 and each nearest neighboring pair is relatively 
rotated in isospin space by vr with respect to the line joining them. At low density, the system 
maintains the original configuration, i.e., it appears as an FCC crystal with the dense centers of 
single skyrmions on each lattice sites. At high density, however, the system undergoes a phase 
transition to a CC crystal made up of half-skyrmions. Here, only half of the baryon number 
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Table 2: Symmetries of the FCC skyrmion crystal 



symmetry 


space 




Uo 




X 


reflection 


X {-x,y,z) 


Uo- 


> (fj, -vri,7r2,7r3) 


Xo 


Xo 


3-fold axis rotation 


X {y,z,x) 


Uo- 


(fj, 7r2,7r3,7ri) 


Xo 


Xo 


4-fold axis rotation 


X {x,z,-y) 


Uo- 


y (fj, 7ri,7r3, -7^2) 


Xo 


Xo 


translation 


X — > {x + L,y + L,z) 




{a, -TTi, -7r2,7r3) 


Xo 


Xo 



carried by the single skyrmion is concentrated at the original FCC sites, while the other half is 
concentrated at the center of the links connecting these points. 

Let us denote the field configuration for the skyrmion field by Uo{x) = a + if ■ tt (with 
a"^ + -k"^ = 1) and Xo(^) for the scalar field. The FCC configuration we are considering has the 
symmetries listed in Tab. 2. There, 2L is the size of the single FCC unit cell that contains 4 
skyrmions. Thus, the baryon number density is p = 1/2L^. Normal nuclear matter density 
occurs at pQ = 0.17/fm^which corresponds to L 1.43 fm. Note that the chiral a field has 
exactly the same symmetries as those of the x- 

As for the constrained fields (cj, tTq), it is convenient to work with ^^unnormalized^ fields 
(fT,7fa), which can then be normalized^ as 

(20) 



and similarly for tTq (a = 1, 2, 3). The field configurations obeying the above symmetries can be 
easily found by expanding the unnormalized fields in terms of Fourier series as |21| 

^(^) ~ l^ahc cos{aTTx / L) COS {biT y / L) cos{cTTz/L), (21) 



and 



TTi{x) = ahki sm{hiTx / L) cos{k-Ky / L) cos{Iitz/L), (22) 

h,k,l 

Ti'2{x) = ahkicos{l7rx/L)sm{h7Ty/L)cos{k'Kz/L), (23) 

h,k,l 

7ts{x) = ^ a/iH cos(/c7rx/L) cos(Z7ry/L) sin(/i7rz/L), (24) 



h,k,l 



and finally 



Xo{x) = 7abc cos(a7rx/L) cos(fe7ry/L) cos{c-kz/L). (25) 

a,b,c 



The symmetries of Tab. 2 restrict the modes appearing in eqs. H21II25P as follows; 

(Ml) if h is even, then k, I are restricted to odd numbers and a, b, c are to even numbers, 

(M2) if h is odd, then k, I are restricted to even numbers and a, b, c are to odd numbers. 

Furthermore, Ohki = Ohik and (3abc = Pbca = Pcab = Pacb = Pcba = (3bac- 

We can locate, without loss of generality, the centers of the skyrmions at the corners of the 
cube and at the centers of the faces by letting a = —1 and 7rj(i = 1, 2, 3) = at those points. For 
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the skyrmion field to have a definite integer baryon numbers per site, we should have a = +1 
and 7rj(z = 1,2,3) = at points such as (L,0, 0). This produces the constraint, 



E ^abc = . (26) 



a,6,c=even 



As for jabc associated with the scalar field, there is no such constraint, but the coefficients should 
be arranged to satisfy Xo ^ 0, a consequence of our choice of vacuum. 

If we had only the modes (IVGJ) in the expansion, the configuration would then have an 
additional symmetry, namely, under the translation x {x + L,y, z) the field undergoes an 
0(4) rotation by vr in the o", vri plane. Furthermore, in order to satisfy the constraint x ^ 0, x 
must vanish identically. This configuration corresponds to the half- skyrmion CC as explained 
above. Because of this additional symmetry, physical quantities such as the local baryon number 
density and the local energy density become completely identical around the points with a = —1 
and the points with a = +1. Thus, one half of the baryon number carried by a single skyrmion 
is concentrated at the sites where the centers of the single skyrmion are expected to be in the 
FCC crystal. The other half of the baryon number is now concentrated on the links connecting 
those points, where the a takes the value +1 and, in the original FCC configuration, the local 
baryon number density is rather low. As a consequence, the expectation value {a) goes to zero, 
signaling the restoration of the spontaneously broken chiral symmetry. Both modes, (Ml) and 
(M2), are included in the actual numerical procedure which we define next. The half-skyrmion 
crystal configuration arises at high density where the expansion coefficients associated with the 
modes (Ml) become small. 

In order to obtain the coefficients we minimize the energy per baryon E/B given by 

E/B = -] f d\CM{UQ,Xo) 

4 J Box 

2 

T,{diUldiUo) + ^Tr \ul^^UQ, uld.Uo 




Tr(2 - Ul - Uo) + -daodixo + V{xo) } (27) 



by taking the coefficients of the expansions as variational parameters. In eq. ()27() . the subscript 
'box' denotes that the integration is over a single FCC box and the factor 1/4 in front appears 
because the box contains baryon number four. We employ "the down-hill simplex method" [H] 
for the minimization process. 

Before going into further details, let us do a rough study of the role the chiral field plays in 
the phase structure of the system. For this, we take xo to be a constant 

Xo/fx = X. (28) 

Then eq.(|27|) can be approximated to 

E/B{X,L) = X\E2/B) + {Ei/B) + X\ErrjB) + {2L^) (x\\nX -\) + \), (29) 

where £"2, E/^ and Em are, respectively, the contributions from the current algebra term, the 
Skyrme term and the pion mass term of the Lagrangian to the energy of the skyrmion system 
and (2L'^) is the volume occupied by a single skyrmion. It can be understood as an effective 



9 




Figure 2: Energy per single skyrmion as a function of the scalar field X for a given L. The 
results are obtained with the (E2/B), {E4/B), and {Em/B) of ref.^ and with the parameter 
sets B. 

potential for the chiral field in medium, modified by the coupling of the scalar to the background 
matter. With the parameter values obtained in ref.,7,] for the Skyrme model without the scalar 
fields, the effective potentials E/B[X) behave as shown in Fig. |21 At low density (larger L), the 
minimum of the effective potential is located slightly from X = \. As the density increases, the 
effective potential V{x) develops another minimum at X = which was an unstable extremum 
of the potential in free space. At L ~ 1 fm, the newly developed minimum can compete with 
the one near X ~ 1. At higher density, the minimum gets shifted to X = where the system 
gets stabilized. 

In Fig. 13 we plot E/B{Xmin,L) as a function of L. The figure in a small box is the 
corresponding value of Xmin as function of L. There we see an explicit manifestation of a first- 
order phase transition. Although the discussion is perhaps a bit too naive, it essentially encodes 
the same physics as in the more rigorous treatment of xo given below. 

Given in Fig. 0] is the energy per baryon, E/B, as a function of the FCC box size parameter 
L. Each point corresponds to the lowest energy crystal configuration for the given value of L. 
The solid circles, solid squares and solid triangles are obtained with the parameter sets A, B, 
and C, respectively. Furthermore, the black solid lines correspond to the single skyrmion FCC 
phase, while the gray lines to the half-skyrmion CC phase. As we squeeze the system from 
L = 6 on, the skyrmion system undergoes at L = Lpt a phase transition from the FCC single 
skyrmion configuration to the CC half-skyrmion configuration. The transition appears to be 
first order. In the half-skyrmion phase, Xoi^ vanishes. The energy of the system comes only 
from the Skyrme term and the scalar field potential. The former roughly scales as ~ 1/L and 
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Figure 3: Energy per single skyrmion as a function of L. The results are obtained by minimizing 
E/B{X, L) with respect to X for a given L. 

the latter exactly as (the volume of the box); explicitly, 

E/B ~ a/L + 2Vixo = 0)L^, (30) 

a being constant. After the phase transition, the energy per baryon E/B continues to de- 
crease (even faster) and reaches its minimum point at L = Lmin- The precise values for the 
phase transition point Lpt and the minimum energy point Lmin depend on the parameters and 
are therefore different for the various sets. 

The incorporation of the scalar field into the Skyrme model Lagrangian produces quite dra- 
matic effects on the properties of skyrmion matter. For example, the energy per baryon drops 
down to ~ 700 MeV. This can be understood since the potential energy between the skyrmions 
amounts to nearly 50% of the skyrmion mass and comes from the medium-range attraction 
generated by the scalar field. 

From here on we denote by (Q) the average value of a quantity Q(f) over the FCC box 
defined by 

(Q) = ^ j^^/rQ{r). (31) 

In Fig. [3 we represent (a) and {xo/ fx) ~ the averaged values of cr(r) and Xo(^) over space - as a 
function of L. Both quantities go to zero at the same phase transition point L = Lpt- However, 
for larger values of m^, the transition properties of the two quantities look different. Note that 
while a{f) ^ locally and (a) = is obtained by the averaging process, (xo/fx) — results 
because Xo(0 = throughout the whole space. We will take (xo) as the effective value of f^, 
i.e., /*. We interpret this result as the vanishing of the "soft" part of the gluon condensate tied 
to the quarkonium component. 
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Figure 4: Energy per single skyrmion as a function of the size parameter L. The sohd circles, 
solid squares, solid triangles are the results obtained with the parameter sets A, B, and C, 
respectively. 




Figure 5: (cr) and as a function of the size parameter L. 
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5 Properties of Pions and Scalar in Skyrmion Matter 

In this section, we study the properties of the pionic and scalar excitations in dense matter. For 
this purpose, we incorporate the fluctuations on top of the static skyrmion crystal as in This 
can be achieved using the Ansatz 

U{x) = ^U^ix)Uo{x)^Ju^{x), x{x) = Xo{x) + x{x), (32) 

where 1/-,^ = exp(z-f • (jj/f-w) as in ref. [7]. Hereafter, we will drop the tilde to denote the x 
fluctuation. 

Expanding the Lagrangian Q up to the second order in the fluctuating fields, we obtain 

^{U, x) = ^m{Uo, Xo) + ^M,x + ^M,4> + ^M,^x^ (33) 

where the subscript 'M' denotes the matter field. Cm the Lagrangian density for the static 
configuration for the background matter and the various terms are given by 



M,. 



M,x<t> 



Here, 



^d^xd'^X - 2 



^y(l + 3ln{xo/fx)]+^4^r{diU^,d,Uo) 



fxJ 



Px 4 



+^l-flml{l-a{x)) 



'X 

,2 



X 



[2X0 fn 

\P 4 



Tr{i{Li-Ri)T'')]xd'(l>a + 



3Xo fTrml 
Px 2 



TT{iTaUo)]xr- 



(34) 

(35) 
(36) 



Gab{x) 



Ht^{x) 



Yj {^ab + 9ab{x)) + 



32e2/2 



1 



2 , 



~b^ 



where for later convenience, we have defined a part of G separately as 

gab{x) = lTr{TaUoTbU^ - TaU) = - {fr'^dab - T^aT^b) 



(37) 
(38) 
(39) 

(40) 



which comes from the current algebra term in the Lagrangian. Here Lj and Ri are defined by 



U = {^^Ul)Uo, R, = {^^Uo)Ul 



(41) 



in terms of the background matter fields. 

The previous equations show how the medium, represented by the background skyrmion field, 
influences the elementary excitations. In order to get an idea of the expected results we show the 
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tree level approximation where we have substituted the background field by its space average 
over the cell. In this way we obtain closed form expressions for the in-medium parameters 
showing their relation to their vacuum values. 

Due to the symmetry of the background field, (Tr(ir(jC/o)) vanishes in the averaging procedure 
and therefore there are no terms with x^Pa in ^M,x<j> ^^^^ order. The relevant terms arise from 
Cm,<i> and Cm,x ^"^^ written as 



72 777*272 
Cm,4> = -fMa-^:^<t>l---, (42) 

1 m*2 

C-M,x = ^XX-^X^---- (43) 

where the pion wave function renormalization constant Z.,^, the in- medium pion mass m* and 
scalar mass m* are defined as 

mfzl = (m^Xa{x)ml\, (45) 




+ ^jlflml{l-a))Y (46) 

The wave function renormalization constant Z-,^ gives the ratio of the in-medium pion decay 
constant /* to the free one, and the above expression arises from the current algebra term with 
/tt in the Lagrangian. The other two equations reflect how the medium affects the effective 
masses of the mesons. 

In Fig. El we show the ratios of the in-medium parameters relative to their free-space values. 
Only the results obtained with the parameter set B are presented. The other parameter sets 
yield similar results except that Lpt takes different values It is important to note that at 
L = Lpt, all the ratios approach zero except for m*/m^. 

Let us analyze these ratios at the light of the above equations and related work. As stated 
before xo and (a) vanish in the dense matter phase. This is the reason for the vanishing of two 
of the above ratios. The non- vanishing of m*/m^ is due to the existence of a pure background 
term which appears in the contribution to the in-medium scalar mass. This term substitutes in 
the dense matter phase the trace anomaly relation eq. by 




^4^i{diUldiUo)), (47) 



which remains finite. In line with our previous arguments, we see that quantities related to the 
condensate in free space become tied to chiral fields in dense matter. 

The extra terms in eq. (|46|) . which are not proportional to , come from the properly scaled 
current algebra and the pion mass terms. In the (matter-free) vacuum these terms describe the 

^Since the parameters are not uniquely given and will strongly depend on the details of the Lagrangian, one 
should not take the precise value of Lpt seriously. 
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couplings of X to the vr's. In dense medium x couples to the static background matter fields 
which contribute to the x mass. 

The vanishing of the pion mass, which already occurred in the pure Skyrme model calculation 
due to the vanishing of the {a) , shows the chiral behavior of xo j since the transition density for 
both phenomena is the same. 

The vanishing of f^/f-K represents the main qualitative difference from our previous calcu- 
lation |7^ and requires a more careful analysis. In our previous calculation the ratio was found 
to saturate to ~ 2/3 at high density. The scalar field shifts the effective potential of the pion 
in matter from one where the symmetry of the ground state is spontaneously broken, to one 
where it is not, i.e., in the language of the sigma model, the scalar field "lifts" in the medium 
the matter-free vacuum constraint . One may interpret the scalar field xo ^ the "radius field" 
of Chanfray et al ^ ^. 

The phenomenon discussed above is closely related to "Brown- Rho" scaling j^. In the de- 
scription of 0, the density dependence comes solely from the change in the mean field x* 
with the corresponding change to the skyrmion structure ignored. Our present result corrects 
and gives a precise meaning to the scaling relation of Similarly eq. (|44|) compares to the 
corresponding f*/ of ref. Again an additional factor, (1 — |7f^(x)), associated with the 
background property appears as a correction. 

At low matter density, the ratio f^/f-K can be fit to a linear function 

^ ~ 1 - 0.24(p/po) + • • • (49) 

At p = pq, this yields f*/fn = 0.76 which is to be compared with f*//-,^ ~ 0.78 of ref. [S]. 
The ratio m^/m^ scales similarly to f^/fn up to p ~ /Oq- In Tab. 3, we list the slopes of the 
ratios. The inset figure in Fig. IHl shows the behavior of the masses m* and m* . They become 
nearly degenerate close to the "critical" density. This near degeneracy may be indicative of the 
"mended symmetry" discussed by Beane and van Kolck ^Hl- Note however that in contrast to 
the dilaton limit of ^B]) the scalar mass remains finite at the phase transition as pointed out 
before while the pion mass vanishes ^. 

It is interesting to look at the decay of x into two pions in the medium. Gathering the terms 
with a scalar field and two pion fields, we get the Lagrangian density for the process x ~^ '^'^ 

^m,xit2 = ^{^ab + gab)xdt^(f>ad''(f)b + ^mla{x)x<pl, (50) 

To compare this result with eqs. (|8I9)) we take only the first term, i.e., the current algebra 
term, in the Lagrangian. Averaging the space dependence of the background field configuration 
modifies the coupling constant by a factor ((xo//x)(l + = ((xo//x)(l ~ |^^))- Taking into 

^It is easy to see this in the work of Beane and van Kolck |16| who relate E oc Ux to the sigma model 
fields a + ir ■ n with of course no sigma model constraint on the radius of the chiral circle. According to this 
interpretation, xo is, trivially, what lifts the vacuum constraint by definition. Furthermore, if we take (xo//x) 
1 + 9/,f,r of ref.^], we can expand /* //^ of Ea. ll44ll as 

/://^~(i+^//^-l^'). (48) 

This is related to eq.(43) of ref. )2U) . (In order for the comparison, we have to take into account that tt corresponds 
to their (f)/ f-K-) Note also that, in ref.|2(J|. the last term proportional to the '^nuclear virtual pion condensation" 
cannot be derived by their "shifted" 8; it is put in by hand. 

®Our skyrmion Lagrangian is not applicable after the phase transition, so the spectra in the Wigner phase 
cannot be considered physically meaningful. 
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Table 3: The slope of the ratios near the origin 





/ttZ/tt 




fx/fx 


ml/m^ 


set A 


- 0.21 


- 0.01 


- 0.14 


- 0.25 


set B 


- 0.25 


- 0.03 


- 0.14 


- 0.28 


set C 


- 0.12 


- 0.005 


- 0.07 


- 0.14 




igure 6: The ratios of the in-medium parameters to the free space parameters. The graph 
small box shows the masses of the pion and the scalar. 
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account the appropriate wave function renormahzation factors, Zjr, and the change in the scalar 
mass, we obtain, in the chiral hmit, an in-medium decay width of the form 



vrvr 



32^ 



((xo//x)(l-ivr2)) 



((xo//x)'(l-ivr2)) 
Applying the naive mean field approximation for averaging, we arrive at 



r*(x 



TTTT 



327r/*2' 



(51) 



(52) 



This expression corresponds to the free one where and are replaced by the in-medium 
quantities m* and /* . We show in Fig. [7| the in-medium decay width predicted with the pa- 
rameter set B. In the region p > ppt where xo = 0) cannot be defined to this order. Near 
the critical point, the scalar becomes an extremely narrow-width excitation, a feature which has 
been discussed in the literature as a signal for chiral restoration |17| I18j . 



6 From an Inhomogeneous Phase to the Half-Skyrmion Crystal 

Up to this point, we have described the background field as an FCC skyrmion crystal. However, 
as can be seen in Fig. 01 in the range of L > Lmin-, the pressure of the skyrmion system 
—dE/dV is negative, which implies that the system is unstable. As discussed in ref. ^2^, the 
system chooses to avoid instability by going into an inhomogeneous phase instead of remaining 
in the homogeneous crystal structure. In this stable phase some part of the volume is occupied 
by the matter with a density Pmin = ^/"^Lmin ^^'^ the rest of the volume is empty. We call this 
"inhomogeneous phase." In this inhomogeneous phase, the spatial average is calculated by 

(g)« = 1 - ((Q)W - {Qtl){p/Pmin), (53) 
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Figure 8: The scaling of the in-medium parameters in the inhomogeneous skyrmion matter. The 
graph in a small box shows the masses of the pion and the scalar. 

where the superscripts '(i)' and '(h)' denote that the average values are evaluated for the in- 
homogeneous and homogeneous phases, respectively, and the subscripts 'vac' and 'min' denote, 
respectively, the vacuum and the lowest energy configuration at p — Pmin • 

What comes out from this naive averaging are the the relevant quantities in the free space 
and those in the minimum energy half-skyrmion phase. Making the averages in eqs. (|44II46() . we 
obtain the simple scaling 

(/v//x)'^*^ = {rjU?^"^ = {mirjUm^f^^ = 1 - p/p^^n, (54) 



m*/m^)2(*) = 1 _ (1 _ (m*/m^)^f2)p/Pmm, (55) 



'X 

and 

where we have exploited that only the ratio ("i* /m;^)^|2 is non- vanishing. 

We show in Fig. |Hlthe schematic scaling behavior of physical quantities in the inhomogeneous 
background matter. Here, the phase transition occurs at p = pmin from the inhomogeneous 
phase to the half-skyrmion phase. The scaling is given simply by ea. (|54l) (which depends only 
on Pmin) and eq. (|55|) . Note that in this naive approximation, the pion mass stays unchanged up 
to 

P — Pmin and then it drops to zero. At low density, we have 

f «f «l-^^WPo) + --- = l-0.18(p/po) + ---, (56) 



(57) 



771 

^«l-0.15(p/po) + ---, (58) 
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where we have used the parameter set B for which Pmin/po ~ 2.8. These results are consistent 
with the scahng found by Song el al ^U]- See eqs. (16-18) of ref. P^ . 

It would be desirable to see how the same quantities scale near the phase transition at 
P = Pmin- Unfortunately, our naive expression does not allow us to do so. 



7 Further Remarks and Conclusion 



This paper represents a significant step toward a systematic understanding of in-medium prop- 
erties of hadrons from the perspective of a unified theory of the hadronic interactions. In our 
previous work ^7 , our starting point was the Skyrme Lagrangian whose parameters are defined 
in order to reproduce mesonic properties in the i? = sector and single nucleon properties in the 
B = \ sector and then we extended the description to a skyrmion matter as a way to understand 
the properties of nuclear matter and the in-medium properties of mesons and nucleons. There 
the phase transition from Goldstone phase to Wigner phase took place from an FCC crystal 
state to a half-skyrmion CC crystal structure. The chiral-symmetry restored phase supported a 
non- vanishing pion decay constant. We interpreted this phase as an analog to the pseudo-gap 
phase of high Tc superconductivity. However there may be a different interpretation. The fact 
that chiral symmetry is restored with a non-vanishing pion decay constant is reminiscent of 
Georgi's vector limit where chiral restoration occurs with the excitation of scalar partners of the 
Goldstone pions with an equal and non- vanishing decay constant [22]. It is known however that 
the Georgi vector limit is not consistent with the chiral Ward identity so we believe that if 
the model is viable, a more likely possibility is the pseudo-gap type realization. 

In this paper, we extended the model studied in j2j to one in which the trace anomaly of 
QCD is implemented in terms of a dilaton scalar field. The idea is very similar to what was 
adopted in [H] and is closely based on the Lagrangian introduced by Ellis and Lanik in 1985 [H]. 

The first question we had to address was how to interpret, in the context of the issue at 
hand, the scalar interpolating field introduced to represent the conformal symmetry breaking 
in QCD in terms of physical particles. It is generally understood that the B = Q sector is 
dominated by pions and gluonium with the trace anomaly more or less saturated by the glueball 
contribution |24l IH] . However chiral symmetry breaking and scale symmetry breaking must be 
intricately tied, requiring the presence of a quarkonium component. How the two components 
are mixed is presently a hotly debated issue and is not our concern here. We simply assume 
that the trace anomaly has two components, the dominant "hard" gluonium component and the 
subdominant "soft" quarkonium component. In a spirit close to that taken by Beane and van 
Kolck [inj, we assume that in matter, the hard component gradually decouples with increasing 
density and in the so-called dilaton limit reached at the critical point, only the quarkonium 
remains active. Thus the "soft" component is instrumental in describing the chiral properties 
in dense systems In order to represent the degree of our ignorance, we consider three possible 
scenarios with the mass of the scalar my- = 550, 720 and 1000 MeV. The results obtained from 
the three do not differ qualitatively, but the densities at which the interesting phenomena take 
place do change substantially. 

The physics we analyze in this new model is closely related with the one presented in our 
earlier work However the details are modified dramatically from the previous work. As 

^In temperature driven chiral restoration, roughly half of the gluon condensate {C^) decondenses across the 
phase transition |25| . Since the gluonium contribution dominates at zero temperature, this means that part of 
the gluonium contribution must also melt at the transition, not just the quarkonium part which is small at T = 0. 
We have no clear idea as to how to unravel this subtlety but we believe this not to be crucial for the qualitative 
argument we are developing here. 
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before, the skyrmion matter we consider has two phases: a low-density phase, which we ulti- 
mately describe by an inhomogeneous phase, and a high-density phase, which is described by 
a CC half-skyrmion crystal. In the former phase which we simulate up to at most a few times 
nuclear matter density by an inhomogeneous structure - the closest we can get at present to a 
Fermi liquid phase - the pions in the medium are massive and their dynamics determined by 
an effective theory where chiral symmetry is spontaneously broken, while the scalar has a large 
mass and couples strongly to the pions. The parameters governing the phase transition, (a) 
and XO) are non- vanishing. In the dense phase, on the contrary, skyrmion matter becomes a 
stable crystal where the two parameters vanish simultaneously at the same density. The pions 
become massless and their dynamics is governed by a chirally restored theory, while the scalar 
remains with a small but finite mass decoupled from the pions. The fact that both parameters 
vanish at the same density confirms the link between chiral symmetry restoration and the "soft" 
component of the x field. The scale anomaly is maintained by the background field and the 
"hard" component, which decouples from the pions and describes scalar gluonium by itself. Our 
model realizes specifically the "lifting" of the radius-field scenario of Chanfray et al j^HI and 
provides a more precise meaning to the scaling behavior proposed in The field xo defines 
the radius of the chiral circle for the in-medium pions and therefore when it is non-vanishing 
their ground state is degenerate and chiral symmetry is spontaneously broken. At the phase 
transition the chiral circle abruptly shrinks to zero and chiral symmetry is restored. In our 
model the background fields, xo and (o", vf), representing skyrmion matter play a dominant role 
in this "lifting" process. 

The model enables one to calculate leading corrections to "Brown-Rho" scaling [5]. The 
in-medium f-,^ decreases with density at approximately the same rate as in jSI for low densities. 
However in our model, two mechanisms participate in this decrease, the changing of xo and the 
deformation of the skyrmion background fields. As mentioned, the main difference with respect 
to our previous work of [7] is that here /^r vanishes in the dense phase whereas it does not in [7] . 
This vanishing is solely associated with the scalar field. 

The masses of the pions and the scalar meson decrease with density in a very characteristic 
way. The pion mass does not change much over the range of density relevant to our consideration. 
Indeed for low density its mass basically remains constant. Not so for the scalar, whose mass 
decreases rapidly, so much so that close to the phase transition it becomes nearly degenerate 
with the pion. This phenomenon may be related to the "mended symmetry" in the dilaton 
limit of the Beane-van Kolck scenario JB] and provides a support to our understanding of the 
behavior of the x field. However the scalar mass remains finite at the phase transition and in 
the dense phase, while the pion mass vanishes identically. The reason for the non-vanishing of 
the scalar mass is associated to the scale anomaly contribution from the background fields. The 
corresponding term is small. However when all the remaining terms vanish according to the 
Beane-van Kolck scenario, the small term becomes relevant. In the true dilaton limit the scalar 
should also become massless and would become the cr-meson of the sigma model Lagrangian. 
However the background obstructs this limit through its contribution to the scale anomaly and 
gives the a a non- vanishing mass. 

We have also calculated the decay width of the scalar in the medium, and we realize that 
at the phase transition the decay width vanishes, signalling the consolidation of the scalar as 
a hadronic stable particle and a possible manifestation of chiral restoration as discussed in the 
literature |T7l[TH| . 

We now make a brief remark on the implications of this work on the parametric dependence of 
an effective field theory matched to QCD discussed in the introduction. What we have computed 
is the density dependence of the ratio ft/ fx and the response to density of the background 
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skyrmion. They should provide density dependence to all variables of interest in the model. 
Since we have no vector degrees of freedom explicitly present in the model, we cannot address 
directly how the gauge coupling g scales. However we can think of implementing the vectors in 
a manner consistent with HLS in which case, their scaling is likely governed by the scaling of the 
the Xo field as first proposed in |Sj with corrections coming from the change in the background 
skyrmion. Thus it is highly likely that we will recover at the chiral transition the limit (7 — > 
together with m*/mp — > as in BR scaling of 

An important feature of the description of the medium properties, discussed in the literature 
|26j - |3nj , is the behavior of the so called pion velocity. This quantity denoted Vj^ is given by the 
ratio of the space part over the time part of the pion decay constant /^//^ and characterizes how 
different effective schemes approach QCD around the chiral phase transition as proposed in |3nj . 
In hot/dense matter the dynamics breaks Lorentz invariance and hence the different components 
could differ in general. For instance, when chiral restoration is induced by temperature, there 
can take place two drastically different phenomena. If one assumes that in the vicinity of chiral 
restoration the only relevant degrees of freedom are pions, then the pion velocity is found ^] to 
vanish — > as T — > Tc. However if there are other light degrees of freedom that can enter in 
the chiral transition as, e.g., in the hidden local symmetry theory of Harada and Yamawaki 
where the light p-mesons with the vector manifestation play a crucial role, the prediction [SOI 
is that the velocity approaches 1 as the critical temperature is approached. There can be a 
small deviation from 1 if one takes into account Lorenz-breaking terms in medium at the scale 
at which the effective field theory (HLS/VM) and QCD correlators are matched. In our present 
scheme where we are concerned with dense matter, the situation is quite analogous to the case 
of HLS/VM. At the leading order of pion fluctuation in the background of the skyrmion matter, 
the pion velocity is 1. This is because the influence of the background at this order is principally 
from the scalar fleld which is Lorentz-scalar. The deviation from 1 can come at the next order 
where the pion interacts with the static soliton background. It turns out that the deviation 
from 1 in v^^ remains small at least in perturbation theory, resembling the case of HLS/VM at 
T ss Tc. The details will be given in a subsequent publication PT] . 

There are several important issues with the model that remain to be addressed. First of all, 
we treated the normal matter to be an inhomogeneous matter to which the FCC crystal state 
collapses. This is not quite the Fermi-liquid state we know nuclear matter should be in. To 
obtain a Fermi-liquid state from skyrmion matter, we need, at the very least, to quantize the 
collective degrees of freedom of the skyrmion matter. This is being worked out at present. The 
solution to this problem will provide us a more realistic information on the intrinsic dependence 
on the background than what we have done here. 

The phase that possesses Wigner symmetry at high density is a half-skyrmion CC crystal and 
it is most likely a solid. It is not clear whether this state - when the color gauge flelds are suitably 
implemented - is color-superconducting as predicted by QCD at asymptotic density |32] or a 
precursor to color superconductivity. Ignoring this, suppose that such a solid state is stabilized 
in compact stars. An interesting question is what the solid state implies for the properties of 
dense stars and their possible observable signatures [HSl • A related question is whether this state 
is different from quark stars or other exotic stars that have been discussed in the literature. This 
is an issue that may be confronted with the phenomenology of compact stellar systems. 
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